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Abstract. The universal analytic expressions in the limit of low temperatures (short 
separations) are obtained for the free energy, entropy and pressure between the two 
parallel plates made of any dielectric. The analytical proof of the Nernst heat theorem 
in the case of dispersion forces acting between dielectrics is provided. This permitted 
us to formulate the stringent thermodynamical requirement that must be satisfied in 
all models used in the Casimir physics. 
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It is common knowledge that dispersion force is a quantum phenomenon which 
results from fluctuating (thermal) electromagnetic fields pp. The most well known 
examples of dispersion interaction are the van der Waals pQ and the Casimir [2-4] 
forces. Dispersion forces play a very important role in surface phenomena, layered 
structures, colloid-substrate interactions, adhesion, foam formation and in stability of 
microelectromechanical systems [5-8]. Recently they were found to be of considerable 
significance in experiments on quantum reflection and Bose-Einstein condensation of 
ultracold atoms near different surfaces [3 Dispersion forces create a free-energy 
difference between materials in the normal and the superconducting phase which may 
influence the value of the critical magnetic field [H]. They are responsible for the 
interaction of atoms and molecules with nanostructures like carbon nanotubes [T2*j . 

The theoretical description of all above-listed phenomena is based on the Lifshitz 
theory [T3|. This theory presents the dispersion force, free energy and entropy between 
the two plates in terms of their dielectric permittivity e(i£) along the entire imaginary 
frequency axis including zero frequency. s(i^) is found by means of the dispersion 
relation using the experimental optical data for the complex refractive index As 
these data are available only within a restricted frequency region, the use of some 
theoretical models of dielectric response becomes unavoidable. Different extrapolations 
of data outside the regions where they are measured (e.g., to low frequencies) may lead, 
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however, to very different theoretical predictions. This places strong emphasis on the 
extrapolation problem in applications of the Lifshitz theory. 

In this paper we present an analytical derivation for the low-temperature (short- 
separation) behavior of the Lifshitz entropy and thermal corrections to the energy and 
pressure between two thick dielectric plates (semispaces) . It is shown to be the same as 
for metals, i.e., universal. We demonstrate that if the plate material at low frequencies 
is described by the static dielectric permittivity, the entropy goes to zero in the limit 
of zero temperature in accordance with the Nernst heat theorem. Alternatively, if one 
includes the nonzero dc conductivity of the dielectric material into the model of the 
dielectric response, the entropy goes to a nonzero positive value when the temperature 
T goes to zero (i.e., the Nernst heat theorem is violated). Finally, we formulate the 
thermodynamic constraint on the extrapolations of the optical data to low frequencies 
and apply it to the topical problem of noncontact atomic friction [15-18]. 

The Lifshitz formula for the free energy of dispersion interaction between two thick 
plates in thermal equilibrium, written in terms of dimensionless variables, is 

^ t) = 32$^ § i 1 - t) r ^ ( x - r " e-y ) + in ( x - r ' e ")i • (i) 

Here, a is the separation between the plates and we use the dimensionless variable 
t = AnksaT '/ '(he) (ks is the Boltzmann constant). The reflection coefficients for the 
two independent polarizations of the electromagnetic field are given by 



etV ~ \/v 2 + Cf(ei - 1) vV + ^ - !) - V 

r -L = / = == • (2) 
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The dimensionless Matsubara frequencies are Q = £j/£ c = It where the dimensional 
ones are £j = 2TTksTl/h and £ c = c/(2a). The dielectric permittivity is computed at 
imaginary Matsubara frequencies E\ = e(i£i) = e(i(i£, c ). 
Applying the Abel-Plana formula [3J 

f (l - i*.) Fit, = f F{t)dt + i f « *%™ , (3) 
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we can rearrange Eq. (pQ) to the form J 7 (a, T) = E(a) + AjF(a, T) where 

he 

32ir 2 a 3 Jo 

/(C, y) = y {in [l - rf (C, y)e~y] + In [l - r 2 ± ((, y)e^] } 

is the energy of dispersion interaction at zero temperature, and 

. w m . iher r°°, F(itr) - F(-itr) . f°°, „, 

is the thermal correction to it. The asymptotic expansions of the energy at both 
short separations and large separations are well known E]- Here we obtain the 
low-temperature (short-separation) behavior of the thermal correction (jHJ) for the case 
of dielectric plates. 
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To solve this problem, it is sufficient to describe the dielectric by its static dielectric 
permittivity Eq = e(0). The reason is that for dielectrics at sufficiently low temperatures 
the Matsubara frequencies giving the leading contribution to Eq. (j^J) belong to the region 
where e practically does not depend on the frequency and is equal to e [this is true 
for AJ 7 but not for E(a)]. To obtain the asymptotic behavior of AjF(a,T) at r C 1 
we, first, expand the function f(x,y), defined in Eq. (JJj), in powers of x = tr. Then we 
introduce the new variable y = y — x to exclude x from the lower integration limit in 
Eq. (jHJ). The subsequent integration of the obtained expansion with respect to y from 
to infinity leads to 



F(ix) 
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where a is real and remains unknown at this stage because all powers of the expansion 
of f(x,y) contribute to its value. Next, we substitute Eq. © in Eq. (jSJ) with the result 



F{a,T) = E{a) 



Tic 
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where C4 = a/240 and ((z) is the Riemann zeta-function. Note that this equation (and 
respective equations for a pressure and entropy) does not allow a limiting transition 
Eq — > 00 in order to obtain the case of ideal metals. The mathematical reason is that in 
our perturbation theory it is impermissible to interchange the limits r — > and 60 — > 00 
in the power expansions of functions depending on e as a parameter. 
The pressure of the dispersion interaction is given by 

dF(a,T) „ , , he 



P(a,T) 



da 
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(8) 
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where Po — —dE/da is the pressure at T = and only the fourth-power term on the 
right-hand side of Eq. (J7J) contributes to the thermal correction. At low temperatures 
this analytical result agrees with the behavior of the Casimir pressure for nondispersive 
dielectrics calculated numerically in Ref. [TT?] . 

Alternatively, the pressure can be found directly from the Lifshitz formula 



P(a,T) 
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Applying the Abel-Plana formula (jHJ in Eq. 0, we get P(a,T) 
where the thermal correction to the pressure is 

iher 



AP(a,T) 
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and the function = $||(x) + $_ 

y 2 dyr\ 
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(x) is defined by 
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By finding the leading term of the expansion of $(x) in powers of x, one arrives at 

.2 



$(ix) - $(-ix) = -i— (y/eE - 1)(£q + £oVto~ - 2) + 0(x b ). 



(12) 
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Substitution of Eq. (fP2|) into Eq. (fTU|) leads to the result 



P{a,T) 
Comparing Eqs. 

Ci 
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and (|T5jl we find the value of so long unknown coefficient 
1 
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Thus, the low-temperature (short-separation) behavior of both the free energy and the 
pressure is given by Eqs. (JJJ), ()13|). (JHJ). By using these results, the asymptotic behavior 
of the entropy of dispersion interaction is described by the expression 

&F(a,T) 3^C(3)(e -l) 2 2 



S(a,T) 



1 



dT ' 64n 3 a 2 (e + 1) 
2vr 2 (5 + l)(e 0v ^ + 2 ^o + 2^ + 2) 



(15) 



135C(3)( V ^+1) 2 

We see from Eq. (jUj) that in the limit r — > (T — > 0) the entropy of both the 
van der Waals and Casimir interactions goes to zero following the same universal law 
which was previously found for ideal and for real metals [Tl]. We have proved that 
the use of the Ninham-Parsegian representation |l j for s(i£i) instead of Eq modifies 
only the terms of order 0(r 5 ) in Eqs. (J7J), (|13|1. The comparison with the results of 
numerical computations for real dielectrics demonstrates that at separations 100-500 nm 
our asymptotic expressions are applicable at T < 60 — 70 K. 

We now turn to a problem of major importance which arises when one includes 
the dc conductivity of the dielectric plates into the model of the dielectric response, 
e~i = Si + Anao/^i = £i + (3{T)/l. Here a is the dc conductivity of the dielectric and 
(3 = Iha^j '{ksT). The conductivity depends on T according to cr ~ exp(— b/T) where 
b is different for different dielectrics. It is significant that for dielectrics the additional 
Drude term is very small for all & ^ 0. For example, f3 ~ 10" 12 for Si0 2 at T = 300 K 
and, thus, it is for sure negligible for all Z > 1. 

One might believe, however, that this term plays a role in the zero-frequency 
contribution in Eq. To test this conjecture we substitute E\ in Eq. (0) and arrive at 



k B T 
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where Li 3 (z) is the polylogarithm function, the asymptotic behavior of T is given by 
Eqs. (J7j), (|14|). and R decreases exponentionally when T — » 0. As a result the entropy 
of the dispersion interaction at T = 0, 
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> 0, 
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in violation of the Nernst heat theorem. Thus, the dc conductivity of a dielectric must 
not be included in the models of dielectric response. This should be compared with the 
case of plates made of real metal (see Refs. [131201 and review |2l a for details), where 
different opinions on the validity of the Nernst heat theorem were proposed. In fact, the 
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mechanisms for the violation of this theorem in some models of metals and dielectrics 
are quite different. In metals, the validity of the Nernst heat theorem depends on the 
scattering processes of free charge carriers on phonons, defects (impurities) etc. For 
the Drude metals with impurities (like in Ref. jTHj) the residual relaxation at T = 
is not equal to zero and the Nernst heat theorem is satisfied. The same takes place 
in the case of metals described by the plasma model. For perfect crystal lattices of 
the Drude metals with no impurities, relaxation at zero temperature is absent and the 
Nernst heat theorem is violated [201 - AH these cases are discussed in Ref. (2Z| devoted 
to metals. For dielectrics, the validity of the Nernst heat theorem does not depend on 
the scattering processes due to quick vanishing of the concentration of carriers when 
the temperature vanishes. Here, the violation occurs due to the inclusion of infinitely 
large dielectric permittivity at zero frequency. Even the sign of the entropy at zero 
temperature for metals and dielectrics is opposite (negative for perfect crystal lattices 
of the Drude metals and positive for dielectrics with included dc conductivity). For a 
complete discussion of this subject, containing all mathematical details, see Ref. [23] 
where Eq. (j!7|) is re-derived in the framework of a more general case of two dissimilar 
dielectrics. 

The above results are important for many applications of dispersion forces. As an 
example we refer to the problem of a noncontact atomic friction where the discrepancy 
between experiment and theory is very large ^3 ^j. In Ref. ^H] it has been 
proposed that the friction observed in the experiment of Ref. [T7j could be due to 
the dc conductivity of an underlying SiC>2 plate described by £/. From the preceding 
discussion, it can be seen that such a proposition would not be in agreement with 
the thermodynamic constraint. Further applications of this constraint in the theory of 
dispersion forces are under way (see Ref. j2H| related to the case of dissimilar dielectrics). 
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